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g : ABSTRACT 



The Bethe-Salpeter equation for pion fluctuations off the chiral soliton in the Nambu- 
Jona-Lasinio model is constructed. By Goldstone's theorem this equation has rotational 
and translational zero modes because the classical soliton is a localized stationary field con- 
figuration which violates rotational and translational invariance. Furthermore, the proper 
^ , normalization of the fluctuating eigen-modes is obtained. Second quantization of the pion 
. fluctuations off the chiral soliton provides an energy functional of the pion fluctuations 
which formally coincides with that of a harmonic oscillator. The corresponding quantum 
corrections to the soliton mass together with the semi-classical cranking prescription yield 
reasonable predictions for the masses of the nucleon and the A-resonance when the con- 
stituent quark mass is chosen to be about 400MeV. These calculations are, to some extend, 
hampered by the non-confining character of the Nambu- Jona-Lasinio model. Comments 
on the 1/Nc counting scheme are added. 
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1. Introduction 



In recent years it has turned out that sohton solutions of mesonic theories provide 
quite a reasonable description of baryons. These approaches are motivated by generalizing 
Quantum Chromo Dynamics (QCD), which is widely believed to represent the theory of 
strong interaction, to an arbitrary number of color degrees of freedom {Nc) In the limit 
Nc oo, QCD reduces to an effective theory of weakly interating mesons (and glueballs); 
the effective coupling constant being of the order 1/Nc ||l|, H- From the fact that baryon 
masses increase linear in Nq, i.e. reciprocal to the effective coupling constant, Witten 
conjectured that baryons emerge as soliton solutions of the effective theory. It is a 
common feature of solitons that their energies increase with the inverse coupling constant. 
These ideas have been brought to practice in the work of Adkins, Nappi and Witten [Q] by 
reviving the Skyrme model 0] from the early sixties. The original Skyrme model consists 
of the non-linear a-model and an antisymmetric forth order (in derivatives) term which 
is mandatory to obtain stable soliton solutions. This model has in turn been applied to 
investigate many properties of baryons ||^. Examples are electromagnetic form factors 
or the phase shift analysis of pion-nucleon scattering [0, Over the past decade the 
Skyrme model then has experienced quite an amount of extensionsQ. E.g. vector mesons 



have been added |]T0[ and the three flavor version of the model has been investigated in 
detaiH. All these extensions were performed such as to comply with the basic symmetries of 
QCD, in particular the chiral symmetry. It has turned out that all these models provide a 
reasonable, not to say successful, description of baryon properties. As a prominent example 
one may quote the natural description of the smallness of the matrix element of the axial 
singlet current, i.e. the famous proton spin puzzle ITT]. Unfortunately, one big problem 



has remained unsolved over the years: The too large prediction for the absolute mass of 
the nucleon. Adopting parameters which are fixed in the meson sector as far as possible 
the nucleon mass is overestimated by about 50%, in three flavor models even more. It has 
only been very recently that quantum corrections to the mass of the Skyrmion have been 
evaluated |jl2|, |13|, |14|. Indeed these have been found to provide a sizable reduction of the 
soliton mass leading to a predicted nucleon mass just of the right magnitude. 

As a matter of fact these corrections, which are due to mesonic loops, can be classified 
by the l/Ai^c-expansion. As the absolute mass of the nucleon (938MeV) is in leading order 
proportional to Nc, the quantum corrections are of one order less, A*"^. Empirically the size 
of these corrections can be estimated from the mass difference between the nucleon and 
the Roper resonance (1440MeV) to be about SOOMeV. The third order in this expansion 
(Nq-^) may be read off from the difference between the nucleon mass and the position of 
the A-resonance (1232MeV). The numerical results in the Skyrme model have been seen 
to follow this pattern. The mass of the Skyrmion is found to be about 1.7-2.0GeV while 
quantum corrections are somewhat less than IGeV |]14|. The mass difference between the 
nucleon and the A-resonance is commonly employed to adjust the only free parameter in 
the baryon sector of the Skyrme model. One may summarize these studies in the Skyrme 
model by stating that the problem of the too large nucleon mass in these models has, to 
a large extend, been solved. What remains to be considered are the analogue quantum 
corrections to other observables. 

The investigations in the Skyrme model have revealed one important fact for the eval- 



uation of such quantum corrections [|T4|: the dominant contribution (about 90%) is due to 



*For a compilation of citations see ref. 
chapter III of ref. g. 
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the zero-modes of the sohton while the contribution of the continuum (or scattering) states 
is almost negligible. These zero-modes arise because the soliton configuration breaks the 
spin (isospin) and translational symmetries of the model. As the quantum corrections are 
naturally subject to subtracting the counterpart associated with the trivial vacuum, which 
does not contain these zero modes but rather possesses a mass gap, it is obvious that the 
zero modes may provide a sizable contribution. The details of this issue will become more 
quantitative in the course of this paper. 

Only a little after the rediscovery of the Skyrme model it has been shown that the 
bosonized versionfl^ of the Nambu and Jona-Lasinio (NJL) model [l^, 0] also contains 



soliton solutions |]18|. Although the NJL model is not as feasible as the Skyrme model the 
NJL model is in a sense superior because it represents a microscopic theory of the quark 
flavor dynamics. The major purpose of the present paper is to examine the corrections 
to the mass of the NJL soliton stemming from quantum fluctuations. In many parts this 



will be similar to the analogous Skyrme model calculations [Tj], however, in the NJL 
model new features will arise since the meson fields here are composite fields built up from 
quarks rather than being elementary as it is the case in the Skyrme model. For example 
the dependence of the action on the derivative of the meson fields with respect to the 
time coordinate, which is essential for the canonical quantization, is not explicitly known. 
As already mentioned it is necessary to subtract the energy of the trivial vacuum (the 
baryon number zero configuration). In the NJL model this provides a further complication 
because in the baryon number zero sector the meson fields do not represent solutions of 
the Klein-Gordon equation as in the Skyrme model but rather of a more involved Bethe- 
Salpeter equation reflecting the quark substructure of mesons. A further motivation for 
the evaluation of the quantum corrections to the mass of the NJL soliton is the fact that 
its classical mass is several hundred MeV lower than the mass of the Skyrmion. It is thus 
interesting to see whether a similar reduction occurs for the corrections and whether this 
will in turn lead to a good prediction for the nucleon mass in the NJL model as well. 

In order to address these questions we have organized the paper as follows. In the 
remainder of this introductory section we describe the NJL model and the appearance of its 
soliton solutions. In section 2 we will introduce small amplitude fluctuations off this soliton 
and briefly review the derivation of the Bethe-Salpeter equation for meson fluctuations off 
the soliton. A thorough study of this equation will provide a suitable normalization of the 
fluctuating modes. In contrast to Skyrme type models this analysis will also be relevant 
for fluctuations in the baryon number zero sector. In section 3 this normalization will then 
prove to be useful in order to derive an energy functional for the meson fluctuations. It will 
be found that these may indeed be quantized using the formalism of second quantization. 
This in turn allows us to apply techniques which have previously been developed in the 
context of the Skyrme model Section 4 is devoted to the discussion of the meson 

fluctuations in the NJL model in channels which contain zero modes. In section 5 we will 
numerically evaluate the contribution of the zero modes to the quantum correction of the 
soliton mass and furthermore make plausible that the contributions of the scattering modes 
are negligible. Concluding remarks are left to section 6. 



As the NJL model is originally formulated in terms of quark degrees of freedom K 
it may be understood as a microscopic model for the quark flavor dynamics. As a matter 
of fact it can be motivated in the large Nc limit of QCD when the gluon propagator 
is assumed to behave like a ^-function in coordinate space []T^|. More importantly the 



NJL model exhibits the chiral symmetry and its spontaneous breaking ||T6|. The latter is 
reflected by a non-zero vacuum expectation value of the quark condensate {qq). In the 
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process of bosonization one introduces composite meson fields which allow one to integrate 
out the quark degrees of freedom [0. The resulting action, ^njl, can then be expressed 
as the sum of a purely mesonic part 

A^= f d^x f-— ^tr(MtM - m°(M + iVft) + {m^f)] (1.1) 
J \ 4Gnjl / 



and a fermion determinant 



Af = Tr log{ip) = Tr log (^i^ - {PrM + PlM^)^ , (1.2) 



I.e. 



ANJL = AF + Am. (1.3) 

Here Pr^l = (1 ± 75)/2 denote the projectors onto right- and left-handed quark fields, 
respectively. As indicated in the above expressions ( p..l| , |1.2|) we will constrains ourselves 
to the investigation of scalar (S) and pseudoscalar (P) fields only. These are parametrized 
by M = S + iP with Sij = S°'t°--/2 and Pij = P"'t°'-/2 representing matrix fields in the 
two dimensional flavor space. Alternatively one may define a polar decomposition of the 
meson fields 

M = (1.4) 

Here the matrix S is Hermitian whereas the matrices C,l,r are unitary. This also gives a 

natural definition of the chiral field U = d^R- ( |1.1| ) Gnjl and m° = diag (m°, m°) 

denote the dimensionful NJL coupling constant and the current quark mass matrix. For 
simplicity we will assume isospin symmetry, i.e. m° = m° =: m°. 

As it stands the action is not well-defined because the fermion determinant diverges 
and hence needs regularization. This is accomplished by first continuing to Euclidean space 
Xq — >■ ir and then considering real and imaginary parts 

' Ar = -Tr ]na f JDr^imJ) 



Ar = -Tr log (pEpE ' j , ^1 = 2^' log i^pEipE ' )-'j . (1.5) 

of the Euclidean fermion determinant separately. Here it is important to remark that the 
Euclidean time, r, has to be considered as a real quantity. Furthermore pE is the Euclidean 
Dirac operator which is obtained from phj the analytic continuation described above. As 
a matter of fact only Ar is divergent and will be regularized by employing Schwinger's 
proper-time description 00]. This introduces a new parameter, the cut-off A, via 



1 r°° ds / +\ 
Ar — ^ --Tr J^^^^ — exp [-spEpE^j . (1.6) 

For a sufficiently large coupling Gnjl the scalar field possesses a non-zero vacuum expecta- 
tion value (S) = diag(m, m) which is related to the current quark mass m° and the quark 
condensate {qq) via the gap equation ||T5 



m = m° + m3^^^^r(-l,(^)^)=m°-2GNJL(gg). (1.7) 
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The quantity m is commonly referred to as the constituent quark mass. Actually the non- 
trivial solutions to eqn ( |1.7[ ) reflect the spontaneous breaking of the chiral symmetry by a 
non-vanishing quark condensate (qq). In order to determine the parameters GnjljA and 
m the fermion determinant is expanded in terms of the meson fields and their derivatives. 
This then yields an effective meson theory which allows one to express physical quantities 
like the pion decay constant /^r in terms of the cut-off A and the constituent quark mass 



m 




f. = ^rol^ . (1.8) 



In practice the physical value /jr = 93MeV is substituted in order to determine A for a 
given constituent quark mass m. Finally Gnjl is determined via Gnjl = m^m/mlf^ where 

= 135MeV denotes the pion mass. Then m is left as the only undetermined quantity. 

In order to examine static soliton configurations it is useful to introduce the Dirac 
Hamiltonian h via 

ipp^ = -d^ - h. (1.9) 

It should be noted that for static configurations only the real part of the fermion determi- 
nant is non-zero. The functional trace in eqn ( |1.6|) is performed in two consecutive steps. 
First one introduces eigenstatesi Qn = (2rz + l)7r/T of dr- Here T denotes the Euclidean 
time interval under consideration. The fermion determinant has been shown |]2T| to acquire 
contributions from the occupation of the valence quark orbit — T£'vai[M] and the polarized 
Dirac sea —TE^^dM]. Whence the total energy due to the fermion determinant is given 
by the sum E^3x[M] + E^i^dM] — E^g^dM = m], which is a functional of the meson fields 
M. Also the subtraction due to the trivial configuration is indicated. Defining eigenstates 
of the Dirac Hamiltonian 

h^, = e,^, (1.10) 



the valence quark part, Eysi[M], may be expressed as ||2T 



E^4M] = NcY.V.\e.\ (1.11) 

V 

wherein t^^, = 0, 1 denote the valence quark occupation number which have to be adjusted 
such that the total baryon number 

^ = E(^--^)sgn(e.). (1.12) 

equals unity. The contribution to the energy due to the polarized Dirac sea is obtained 
from Ar in the limit T — oo. Then the temporal part of the functional trace reduces 
to a spectral integral while the remainder of this trace is performed by summing over the 
eigenstates of h. These manipulations result in [^] 



'^^n are the Matsubara frequencies. The fermionic character of the quarks requires anti-periodic bound- 
ary conditions. 



5 



The classical soliton energy is finally given by the sum 



Eei[M] = E^M] + E,^[M] - E.ac[M = m] + E^[M] (1.14) 

wherein Em[M] originates from the mesonic part of the action ( |1.1D . 

To be specific we employ the hedgehog ansatz for the chiral field (in unitary gauge 

eo(r) = exp Qr • f e(r)) (1.15) 

while the scalar fields are constrained to the chiral circle, S = (S) = m. Substituting this 
ansatz into the Dirac Hamiltonian ( |1.9|) yields 

ho = a. ■ p + [3m (cos9(r) + i'^^T ■ r sin6(r)) . (1-16) 

This Hamiltonian possesses the celebrated feature to commute with both parity (fl) and 
grand spin operators (G = I + cr /2 + t /2). The latter is the sum of orbital angular 
momentum (/), spin (cr/2) and isopin (t/2) operators. Thus the eigenstates '^i, fall into 
independent subspaces which are characterized by the eigenvalues of 11 and G. For the 
hedgehog ansatz the mesonic part of the energy is given by 

Err^Q] = ^T^mlfl I drr^ (1 - cos0(r)) . (1.17) 

The stationary condition 6Eci[Q]/6Q = is obtained using the chain rule. First the 
energy functional is differentiated with respect to the one-particle energies Subsequently 
these eigenvalues are functionally differentiated with respect to 0(r). This leads to the 
equation of motion 

dQ p^{r,r)i-^5T -r = sine(r) hi J dQ {r,r) - — ^| (1.18) 

where the traces are over flavor and Dirac indices only. According to the sum ( p..l4p the 
scalar quark density matrix p^{x, y)) = {q{x)q{y)) is decomposed into valence quark and 
Dirac sea parts: 

p%x,y) = p^^^{x,y)+ p^^^{x,y) 



r vac \ 



^^y) = ^E^-Nerfc(|^|)^'.(2/)sgn(e,). (1.19) 



For various boundary conditions the explicit form of the eigen-functions "^ui^) may e.g. 
be found in refs. p^ , p^ , p^ . Numerically the soliton is obtained iteratively. One starts 



off with a test proflle G(r) which is employed to compute the eigenvalues and -states of hg 
( |1.16| ). These are then substituted into the equation of motion ( |1.18| ) in order to update 
the proflle function B(r). This procedure is repeated until convergence is reached. 



^Here we attach the subscript "0" for later reference when fluctuations are included. 
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2. Mesonic Fluctuations off the Soliton 



Mesonic fluctuations off the chiral soliton in the NJL model have been examined at 
length in the context of the bound state approach to the description of hyperons . 
As the general formalism is not altered we will only point out the key steps and emphasize 
the changes due to the two flavor reduction. 

Here we will only consider pseudoscalar fields and constrain the scalar fields to the 
chiral circle. Then a convenient parametrization for the fluctuations off the soliton is given 
by H 



M(r, t)=m Ur)^f{r, t)^f{r, t)Ur) (2.1) 

wherein ^o(''") denotes static hedgehog configuartion (|1.15| ) while S,f{r, t) contains the small 
amplitude fluctuations r}{r,t) 

0(r, t) = exp Qr7(r, t) ■ r) = 1 + '-r]ir, t) ■ r - ^-qir, t) ■ r]{r, t) + . . . . (2.2) 

Subsequently the total action is expanded up to quadratic order in r}{r,t). The zeroth 
order term just renders the static energy function which is explained in the introduction 
while the linear term vanishes subject to the stationary condition ( |1.18| ). The expression 



for the quadratic part has been derived in ref.[^. In order to present that result it is 
useful to define the perturbation of the Dirac Hamiltonian 

i(3pE = -dr-h = -dr -{ho + hi + h2 + ...) (2.3) 

wherein the subscript labels the power of the meson fluctuations. The zeroth order part 
is already presented in eqn ( |1.16D . In order to display the linear and quadratic parts we 



make use of the unitary matrix |15, 28, 26 



T = ^oPl + ^IPr (2.4) 
which contains the information about the static soliton. Then 

hi{r, —it) = imT/S'y^r] ■ tT' and h2{r, —it) = ——T(3ri ■ tiTK (2.5) 

Here it should be remarked that the temporal argument of the fluctuations is continued to 
Euclidean space, i.e. rj = r]{r, —it) The defining equation ( |2.3| ) is then employed to expand 

the operators J/^ePe and pE ( -pl ) up to quadratic order in r). The resulting expression 



is subsequently substituted into eqns ( |1.6D and ( |1.5D for computing the contribution of the 
fermion determinant to the part of the action being quadratic in the fluctuations. At this 
point it is important to note that one is forced to start off the expansion at eqn (p,.6D rather 
than employing standard perturbation techniques to the energy eigenvalues (|1.10|) . The 
latter prescription would give incorrect results since [dr, /ii,2] 7^ 0. For the expansion in 
Euclidean space one Fourier-transforms the fluctuations 

r7,(r, -it) = ^Ur, ioje)^-'^^^ ■ (2.6) 

J-oo /TT 
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This transformation may be transferred to the Hamiltonian: 

/■+°° dujE ~ _ ■ 

hi{r,-iT) = / ——h^i){r,iujE)e ^""-^^ and 

ft.(r,-.r) = /^~^/^°°^ft,„(..,.,.,.ye-'<--i)' (27) 

because T is static. The final expression for the harmonic part, A2, of the Euchdean action 
is then continued back to Minkowski space: iuje ^- The result may be expressed as a 
sum of three pieces 

A2 = ^val + Aac + Am- (2.8) 

The contribution due to the exphcit occupation of the valence quark orbit (|val)) is given 
by 

Aal = -r/™'iVc / — (val|/i2(r,cu,-a;)|val) (2.9) 

+ H (val|/ii(r,a;)|/i)(^|/ii(r, -^)|val)- ^^'"^ 

/.^vai l^vai - e^j -uj 

The contribution from the polarized Dirac sea is somewhat lengthy 
Nc [°° ds _„2 /■+~ duj 



Aac = — ^=J22e^e ''^ —{fi\h2ir,uj,-uj)\fi) 
2 J1/A2 ^/ATCS ^ J-oo Ztt 

Nc , rs~,^ /■+°° duj 



J1/A2 V47r^J-oo 27r 

" ^ + [uJ^ - (e^ + e,f]R{s- cu, e^, e,) i . (2.10) 

The regulator function R involves a Feynman parameter integral which reflects the quark 
loop in the presence of the soliton 

R{s;uj, e^, eu) = J dx exp (^—s[{l — x)e^^ + xel — x{l — x)uj'^]^ . (2-11) 

Here it should be noted that, upon explicit calculation, it can be shown that the matrix 
elements {^\h2{r, u, —u)]^) and {^\hi{r, uj)\h'){h'\hi{r, —uj)\^) are invariant under u — >• —uj. 
This is the reason why the terms of odd powers in u have been dropped in the eqns 
(|]|,|1|). Stated otherwise: The imaginary part of the action does not contribute in 
the two flavor reduction contrary to the case when strange degrees of freedom are present 
|2^ . Finally the action is completed by the mesonic part 

Am = I'^ ^'l ^ ^ f){r,uj)-fi{r,-uj). (2.12) 

Formally the harmonic part of the action can be expressed with the help of local ^"^(r) 
') kernels 

A2 = y^^^l^ld'rJd'r'^t{r,r',co)Ur,uj)v,ir',-u) 

+ Jd''r<l>l\r)Ur,uj)fjb{r,-u)\. (2.13) 



and bilocal ^f'{r,r',uj) kernels 



These kernels can be extracted from the above expressions. The local kernel turns out to 
be diagonal in isospace $"^(r) = $i(r)5"^. Here we recognize the usefulness of the above 
defined unitary matrix T ( p. 41 ) because it considerably simplifies the presentation of these 

kernels by defining chirally rotated wave-functions \E'^. Then we find 

= -mlf^cosQ{r) + 2r]^^iNcm^l^^{r)f3^,^i{r) 

ds 



2Ncm r^^Yl e,e-l^l{r)(3^,{r) (2.14) 



'1/A2 y/A-KS 

and 



Nc 2 
m 



^",^^\/?E^I(^)/?75r'^^^(r)^t(^')^^5/^^(^') 
1/A2 V 47r 

X + [iv^ - (e^ + e,)2]i?(s; u, e^, e,) . (2.15) 



In ref. it has been demonstrated that these kernels are diagonal in parity and grand 
spin, i.e. meson fluctuations with different grand spin and/or parity quantum numbers 
decouple. This is a consequence of the fact that the classical soliton carries zero grand spin 
and has definite parity. This decoupling of the meson fluctuations will later on be helpful 
since it allows us to consider rotational and translational zero modes separately. The 
equation of motion for the fluctuations is finally obtained by varying (|2.13| ) with respect 
to fja{r,uj) 

I £r' $f (r, r', uj)fi,{r\ u) + $i(r)r/,(r, u) = (2.16) 

which in fact represents the Bethe-Salpeter equation for the pion fluctuations in the soliton 
background. This equation has the following interpretation: The frequency uj has to be 
adjusted to uji such that ( |2.16| ) is satisfied for a non-trivial fjair.oji). The frequency uji is 



called eigen-frequency and fia{r, uoi) denotes the associated eigen-mode or -wave-function. 
Below we will explain how the boundary conditions, which are imposed on the Dirac 
spinors ^'y, transfer to the meson fluctuations and subsequently lead to a discrete meson 
spectrum as well. The numerical methods which are used to solve eqn ( ^.16| ) are reported 
in refs.m, 



Later on we will have to evaluate overlap matrix elements between the solutions to the 



Bethe-Salpeter eqn ( p.l6| ) and states which solve this equation in the absence of the soliton 
i.e. 9 = 0. The latter states are obtained by computing the above defined kernels $i and 
$2 as mode sums which contain the eigenvalues and -states of the free Dirac Hamiltonian 



Let us denote the resulting kernels by (^'^!^'' and the solutions to the corresponding Bethe- 



hi,ee = oc ■ p + (3m. (2.17) 

is by ^^''2^^ and the solutions to the corres 
Salpeter equation by fi^'^\r,ujf'^). Of course, the eigen-frequencies uf^^ will differ from 
those obtained in the presence of the soliton. It should be remarked that 
not necessarily solve a Klein-Gordon equation due to the composite nature of our meson 
fields. 
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In the first place we therefore have to find the proper normahzation of the solutions to 
eqn ( p.l6| ). In the case of the NJL model this turns out to be more involved than e.g. for 
the Skyrme model because the the frequency uj appears with all even powers in eqn ( |2.16|) 
rather than only quadratically. Therefore the metric which appears in the scalar product 
of the underlying Hilbert space of the fluctuations unavoidably depends on the frequency. 
In order to compute this metric we (formally) expand the Bethe-Salpeter equation ( |2.16| ) 
in terms of the frequency 



oo „ 

/dV0f(r,r')a;^r^,(r',c.,) = O 



(2.18) 



and assume that fja{r,u!i) represents an eigen-mode with frequency Ui. The expansion 
( p. 18 ) obviously represents a generalization to the free Klein-Gordon equation which cor- 



responds to Of) 



-ml 



d and Oi = 1 while all other 0„ vanish. Also in the Skyrme 



model only Oq and Oi are non-zero. As compared to the Klein-Gordon equation, however, 
they acquire additional space dependent factors. 
For the general case we start with the identity 







(2.19) 



Noting that the O!^ are Hermitian under the spatial integration we obtain by substituting 
the Bethe-Salpeter equation ( |2.18| ) 



0= d\ d\'J^Ur,uJ.)0:^'{r,r') 



n=l 



2n 2n 



{u^-u^) f dh f d'r'Ur,u;,) 



n-l 



2p 2{n-l-p) 



n=l 



p=0 



(2.20) 



f2.211 



Here the relation — b'^ = {a — b) J2m=o o.'^b^"' ^ has been used. Eqn ( |2.21| ) allows one 
to impose the orthonormalization condition 



d'r / dVf]air,uJi) 



n-l 



p=0 



2p 2(n-l-p) 



n=l 



r]b{r',Ljj) = 6ij. 



(2.22) 



In particular a solution ria{r,uji) to the Bethe-Salpeter equation is normalized to 



d'r / d\'Ur,^^)M''\r,r',u,)fib{r',iv,) = 1 



which defines the metric tensor 

d^2^{r, r', to 



n=l p=0 
ab I 



n— 1 oo 
2p 2(n-l-p) 



n=l 



2{n-X)y^ah 



Ol\r, r') 



(2.23) 



(2.24) 



As a matter of fact this metric tensor can easily be obtained from eqn ( |2.15| ) in contrast 
to the expansion coefficients O^*. We have thus succeeded in deriving a normalization 
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condition for the solutions to the Bethe-Salpeter equation in the sohton background. This 
represents a special achievement because in the Bethe-Salpeter equation the frequency u) 
appears with arbitrary even powers. In coordinate space these are time derivatives. 

As we will take advantage of the grand spin symmetry when computing overlaps be- 
tween r) and fj^^'' it is appropriate to also define a metric tensor }Af'^^\r^r\ujf'^) for the 
fluctuations in the baryon number zero sector. This quantity is obtained by substituting 
$f in eqn (pf. 

Now we are finally enabled to equip the overlap {f] (r,u;i)|77(°)(r,cuf )) with an ingenious 

meaning. To do so we first redefine the meson wave-functions f] ^ (p such that they are 
normalized to unity with respect to the trivial metric 

J £r(f>{r,uj,)-cl>{r,cui) = l. (2.25) 

This can obviously be achieved by the introduction of the root |^ (see section 4 for its 
exphcit construction) 



M''\r,r',u,) = J d'xY.(VMy\r,x,Ui) (v^)'' (r', ^, u;,) (2.26) 

c=l 

into the wave-function 

0»(r,^,) = J d''xY,[VMy\x,r,uj,) Ux,u,). (2.27) 



c=l 



The application of these definitions to the baryon number zero sector is straightforward 
yielding (f)^^\r,uj^^'^). To this end the above mentioned matrix element is defined as 

{fl{r,Ui)\fi^'\r,uf)) := J d\ cP{r,u,) ■ (P^'\r,u;f). (2.28) 



m 



This actually is the generalization of the Skyrme model definition for this overlap 
the case that the metric in the orthogonality condition ( p.22| ) depends on the frequencies of 
the states. It should, however, be noted that this procedure only yields normalized wave- 
functions and that there does not exist an orthogonality condition for the fiuctuations 
(f){r,uJi). Nevertheless the definition ( p.28| ) represents the most reasonable one because the 
wave-functions ff and f)*-"-* obey Bethe-Salpeter equations which in principle are discon- 
nected because they belong to different baryon sectors and thus different Hilbert spaces. 

There is one more important point which has to be mentioned in the context of the 
Bethe-Salpeter equation for the meson fluctuations in the NJL model. As this stems from 
a shortcoming of the NJL model in general it effects the eigen-modes in the presence as 
well as in absence of the soliton. The NJL model is well known not to possess quark 
conflnement. Thus meson flelds may decay into "free" quark-antiquark pairs once the 
meson energy ui is beyond the two quark threshold (^thLs ~ ^m. Technically this appears 
because the argument of the exponent in Feynman parameter integral ( |2.11| ) turns neg- 
ative. As a matter of fact the analytic continuation from Euclidean to Minkowski space 
yielding eqn (|2.1(]| ) is ill-deflned for > wlhreg. Along this path in the complex plane the 
logarithm develops an imaginary part which measures the width for the meson to decay 
into a quark-antiquark pair. Therefore the expansion of $2'''''''* which flnally yielded A4°''''°'' 
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converges only for u^^^ < 2m. Later on we will have to perform sums over the eigen-modes 



r/'^°^(r, cuj^"*). As a consequence of these considerations, euchres provides a natural cut-off to 
these mode sums. In the presence of the soliton the situation is even worse because the 
valence quark orbit gets bound and acquires an energy eigenvalue evai < m. Following the 
analysis presented in appendix B of ref.[^ this leads to the threshold cuthres = 2|evai| < 2m. 
Fortunately, this threshold will not be of utmost relevance for the ongoing studies because 
we are mostly interested in the zero modes (wj = 0) in the soliton background. 

3. Energy Functional for Fluctuations 

In order to compute the quantum corrections to the soliton mass we first have to con- 
struct the energy functional for the meson fluctuations. Again this is not straightforward 
because the Bethe-Salpeter equation involves arbitrary even powers of the frequency uj. We 
therefore have to go back to the expansion ( p. 18 ) of the Bethe-Salpeter equation. Formally 
the Bethe-Salpeter equation is obtained from the action formulated in Fourier space 



(0) 



.n=0 



S[ff] = lj^Jd\Jd'r'Ur,-. 

Undoing the Fourier transformation one obtains the Lagrange function 

1 /■ r 



(3.1) 



(3.2) 



n=0 



where the superscript denotes the order of the time derivative of the fluctuation ria{r,t). 
As in classical mechanics we compute the total derivative of L with respect to the time 
coordinate in order to derive the conserved energy. Upon integration by parts we find 



d 



E 

n=0 



d^r 



' d 
It 



'II 



6L 



m=0 



+ (-l)-r^W(r,t) 



5L 1 



(3.3) 



Inserting the expansion 
last, i.e. surface term in 



2]) and the Fourier-transformation for the fluctuations (|2.6|) the 
]3D is shown to vanish for meson fields which satisfy the Bethe- 
Salpeter equation (|2.16| , p.l^ ). Thus the conserved quantity which has to be identified as 
the energy functional and thus the Hamiltonian is given by 



d^r 



1 n=0 m=0 



6L 



5r,'t\r,t). 



„ „ ( oo n 

J d'r J d\' Y: Ei-^rV^^'"'\r,t)0:,\r,r')4-'-"'\r',t) 



n=0 m=0 



o E^i"Hr',t)Of(r,r')r/(")(r',t) . 



(3.4) 



m=0 



Next the fluctuating field rf{r,t) is decomposed in terms of the solutions to the Bethe- 
Salpeter equation ( p.l6| ) 



(3.5) 
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because the solutions to ( |2.16| ) come in pairs icUj, i.e. r]{r,t) describes a real field. The 
canonical quantization prescription then corresponds to require the commutation relations 



J 



6ij. 



(3.6) 



The decomposition (|3.5|) is substituted into the energy functional for the meson fluctuations 



The aim is then to obtain an expression in terms of the operators and a- by using 
the orthonormalization condition ( |2.22| ). It is obvious that at an intermediate step the 
term involving has to be eliminated. This is done in a symmetric way by the help of 
the Bethe-Salpeter equation 



^ ■' ■' Ln=l 



• (3.7) 



After a somewhat tedious calculation we find 



2{n-l)y^ab 



+E 



X 



oo 
.n=l 

oo n— 1 



d^r / rfV r]a{r,uji) 



.n=l m=0 



fib{r',ujj) 



\2 i(uji+u)j)t 



. (3.8) 



Here we finally recognize the appearance of the orthonormalization condition (|2.22|) . Hence 
we arrive at the energy operator of an harmonic oscillator 



= E^i ( + 



(3.9) 



Although this result is not unexpected it is at the same time non-trivial because the Bethe- 
Salpeter equation (|2.16|) , which is the defining equation for the normal modes of energy 
uji, contains arbitrary even powers of the time derivative operator when transformed to 
coordinate space. Needless to mention that the above analysis goes through as well in the 
absence of the soliton resulting in 



= E 



, ,(0) /^^(o)t.(o) , 1 



(3.10) 



Here af^^ and af^ respectively denote the creation and annihilation operators for the 
eigen-modes f]^'^\r .,ujf^). 

We have thus seen that the vacuum energy, Y.i^il'^i corresponds to the one of an 
harmonic oscillator and that this result is independent of the specific form of the kernel 
for the Bethe-Salpeter equation. There are two restrictions only: 
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(i) The eigenvalues appear in pairs icjj. 

(ii) The background field is static. 

In order to see that the vacuum contribution to the Hamiltonian ( p.9| ) comes with no 
surprise but rather is a feature of the semi-classical treatment of the meson fiuctuations 
we read off the inverse propagator for these fluctuations from the action ( |3.1D 

, oo 
ri=0 

Since the background fleld is static the propagator V is local rather than bilocal in the 
frequency u. 

In order to extract the vacuum contribution of the meson fluctuations in the soliton 
background in the semi-classical approximation one considers the functional integral (con- 
tinued to Euclidean space) 

(3.12) 
i.e. 

= ^ Tr log (P-I) . (3.13) 

The temporal part of the functional trace can be done because the propagator is local in 
the frequency. As in the case for the fermion determinant (c/. section 2) the vacuum part is 
obtained by considering the limit of inflnitely large Euclidean times T. Then the temporal 
part of the trace becomes a spectral integral 

^M = |/_^^ Tr\og{v-\u:')) (3.14) 

where we have suppressed spatial and isospin arguments. Furthermore Tr refers to the 
trace over all degrees of freedom other than the time coordinate. Next we integrate by 
parts 

The surface term has disappeared since the propagator only depends on uj"^. The trace Tr 
can be evaluated with the help of the solutions to the Bethe-Salpeter equation f){r,uji) 

Am = \ jy^2u^Y.^mu')^^V-\u')\^. (3.16) 

Now it is important to note that the states \i) represent the solutions to the Bethe-Salpeter 
equation. This allows one to expand 

{i\V{uj') = -^A^,m-\u:')\ +0{u:'-ujt) \ \ (3.17) 



Dt] exp{-5'£;[r7]} 



Det I? 



1-1 
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Hence 



T f du ^ 2uj'^ 
Am = -77 / 7^ 



2n - ujf 



0(u;^--u;!)\ ' l-V-\„^)\i)^ (3,18) 



Under the assumption that summation and integration may be exchanged the spectral 
integral may be computed with the help of the residue theorem 



= -77 E 1™ (^-^.)^ -^{i\\j—V-\u') V'\u:^)\i) 



' d 



= -|E^- (3-19) 

Thus we have seen on the formal level that the above listed restrictions to the Bethe- 
Salpeter equation are sufficient to provide a vacuum energy which is given by a sum of 
eigenfrequencies. However, the above derivation has to be taken with some care because 
the functional traces are, of course, divergent and thus require regularization as do the 
expressions ( p.9|) and ( p.lO|) . This issue will be discussed next. 

The main content of the above considerations is indeed the fact that we have obtained 
a quantized Hamiltonian for the meson fluctuations which are formally equivalent to a 
harmonic oscillator. Considering eqn ( |3.19| ) one might be tempted to identify the quantum 
corrections to the energy as the difference 

2 ^ 2 ^ 

However, this expression diverges logarithmically. In the context of the Skyrme model 
Holzwarth has shown that two additional subtractions yield a flnite (renormalized) result 
[|T^. To identify these subtractions in the NJL model we make use of the fact that eqns 
( p.9| ) and ( p.lO[ ) imply the existence of operators and Hq which have the eigenvalues 
Ljf and oof'^'^ when acting on the meson eigen-modes in the baryon number one and zero 
sector, respectively. One may especially deflne the "perturbation potential" V via 

= Hi + V. (3.20) 

A finite expression for the vacuum energy is obtained after subtracting the first three terms 
in the expansion 



Tr [^jHl + V) = Tr {H,) + -Tr [H,'V) - -Tr [H,'VH,'V) + . . . (3.21) 
resulting in the finite energy correction due to the quantum fiuctuations |I^ 

^E=]^Tr[H-H,- ii/o" V + ^^0"'^') • (3-22) 

This expression, which was derived by Holzwarth in the case of the Skyrme model, repre- 
sents the actual starting point of our investigations in the context of the NJL model. The 
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main complication in comparison with the Skyrme model is that we do not have explicit 
access to the operators H and Hq. However, we are able to compute the associated energy 
eigenvalues by solving the Bethe-Salpeter equation (p.l6|). Furthermore we have a suitable 



definition ( |2.28| ) for the overlap of the eigenstates of H and Hq. The energy correction 
( |3.22D actually represents the 3 + 1 dimensional generalization of the result obtained by 
Cahill, Comtet and Glauber [P2| in a 1 + 1 dimensional model. It should, however, be 
noted that the correction ( ^.22 ) is not free of ordering ambiguities. 

As already mentioned we have access to the eigenvalues of H and Hq only. Then the 
correction ( p.22| ) may be expressed as 



i I .7 



(3.23) 



X 









4" 




3 + 6 








1 


a- 









r.uy)) shows up since the trace has to be performed in a 



Here the overlap {fj{r,uj. 
common Hilbert space to chop off both frequency sums in the same way. At this point it 
should also have become clear that the formal considerations in section 2 on the orthonor- 
malization of the solutions to the Bethe-Salpeter equation are in fact unavoidable. Eqn 



( |3.23|) obviously yields AE = in the absence of the soliton, i.e. when Ui 



UJ. 



(0) 



and 



(0) 



). We furthermore observe from eqn ( |3.23|) that AE is of the order 



A*"^ in agreement with the assertions made in the introduction. 

From eqn ( |3.23| ) it is obvious that the zero modes, i.e. states with cuj = may lead to 
a sizable reduction of the total energy E + AE since for these states AE is negative. The 
fact that for scattering states Ui ^ uif'^ demonstrates that the main contributions to the 
quantum corrections is indeed due to the existence of zero modes in the soliton background. 
This can also be understood in the context of the phase-shift 6{k) expression for the Casimir 
energy, E^as ~ ("V^^r) / dk6{k) up to counterterms which render this integral finite for 
k ^ oo [Q. According to Levison's theorem an additional vr has to be included for each 
bound state. The only bound states in the background of the NJL soliton are the zero 
modes. Thus the channels in which these modes appear may considerably contribute to 
i^cas- In the proceeding section we will therefore concentrate on the channels which contain 
the zero modes in the NJL model. 

4- Zero Mode Channels 



Zero modes, i.e. solutions to the Bethe-Salpeter equation ( ^.16| ) with uji = 0, arise 
whenever the stationary background field (the "vacuum" seen by the meson fiuctuations) 
breaks a symmetry of the underlying theory. In that sense the zero modes are Goldstone 
bosons. In the case of the chiral soliton the hedgehog field configuration (|1.15| ) violates the 
rotational and translational invariance. We therefore expect zero modes to be associated 
with infinitesimal spatial rotations and translations of the soliton. Due to the grand spin 
symmetry the zero mode corresponding to infinitesimal iso-rotations is equivalent to the 
one of the spatial rotations. Although the model is invariant under axial rotations (for 
massless pions) a corresponding zero mode does not exist since the infinitesimal axial 
rotation does not leave the vacuum configuration (M = m) invariant. 
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In order to identify the zero modes we expand the parametrization ( p.l| ) up to hnear 
order in r}{r, t) 

M = m {^o' + i^oV ■ r^o} + . • • = Mq + im^ot] ■ r^o + • • • • (4.1) 

For the extraction of the formal structure of the zero modes we have to identify the hnear 
term with [Q, Mq]. Here Q refers to the generator of the symmetry transformation. For the 
spatial rotation this gives 

r]ji{r) = sin6(r)f x (4.2) 
where 5r is a measure for the infinitesimal rotation. Similarly the translation defines 

rj^(r) = e'(r)rr ■ 5t + {5^ - rr ■ 5t) ■ (4.3) 

r 

It is straightforward to verify that both 77^ and rjrp carry unit grand spin, i.e. these are 
dipoles in grand spin space. In accordance to the Skyrme model notation ||^ we will refer 
to the channel which contains the rotational zero mode as magnetic dipole (Ml) while the 
channel with the translational zero mode is called electric dipole {El). 

Prom the consideration of the zero modes we have obtained the quantum numbers of the 
fluctuations in the Ml and El channels. This allows us to make ansdtze, which separate 
the radial and angular dependencies, for general fluctuations in these channels 

■ iMiir, uj) = T-{rx C(r, u)) = ^[r ■ C{r, w), r ■ f] (4.4) 

and 

T ■ flEiir, u) = T- CA^r, uj) + T-rT ■ Csi^, uj)t ■ r. (4.5) 

Altogether these ansdtze introduce nine radial functions. Although the parametrization 
( [4.5| ) is not intuitively clear from ( [4.3| ) it is the most convenient one since the action of the 



grand spin zero object 75T ■ r on the quark wave-functions "^ui^) is well-known and 
does neither change grand spin nor parity quantum numbers. This knowledge is also the 
reason why we expressed ffMi terms of a commutator. In order to compute the matrix 
elements of the ansdtze ( |4.4| ) and ( [4. 51 ) we thus only require the matrix elements of r times 
a radial function. 

In the above notation the zero modes are parametrized by 

C-(r) = sine(r) Sr (4.6) 

Cr-M = l(mr) + '-^)sr. CrW = i(e'W-2^!^)*.. (4.7) 

Now the main task is to substitute the ansdtze ( [4.4|) and (|4.5| ) into the action functional. 
Although these calculations are straightforward they are quite tedious and we do not go 
into the details here. As a matter of fact the presentation of the associated formulae 
would approximately double the length of this paper. Let us rather display the generic 
form of the Bethe-Salpeter equations for radial functions deflned above. As a matter of 
isospin invariance the action only depends on the combinations (^{r,uj) ■ Ci^'j for the 
Ml channel on the one hand and <^^(r, cu) ■ ^^(r', — cu), Cj^{r,u!) ■ ^^(r', — cu) as well as 
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(^g{r,uj) ■ C_b(^'5 foi' the El channel on the other hand. Thus the isospin invariance 
reduces the number of independent radial functions from nine to three. We label these 
by ({r,u), (A{r,u) and ^^(r, u;) and keep in mind that they are three-fold degenerate. 
This causes an overall factor 3 in eqn (|3.23|) when one considers the Ml and El channels. 



The Bethe-Salpeter equation for C('")'^) becomes an homogeneous integral equation in the 
radial coordinate 



/ rfrV2$fi(r,r',cu2)C(r',cj) + $fi(r)C(r,cj) 







(4. 



while the Bethe-Salpeter equation for radial function in the electric channel couples CA{f, t^) 
and CB{r,uj) 




dr r 



+$fi(r 



(4.9) 




dr r 



+$f (r 



CB{r,uj) - -CA{r,uj) 



(4.10) 



Upon explicit computation it can be shown that the non-diagonal elements of the Bethe- 



Salpeter kernel satisfy $ 



ElAB, 



r, r 



ElBAi 



^) which, of course, reflects the 



Hermitian character of the Bethe-Salpeter kernel. As a further consequence the diagonal 
elements as well as ^^^^^ and turn out to be symmetric. 

The numerical treatment of equations like 



is described at length in ref.|30]. Thus 
we will only explain the key steps. As the diagonalization ( [1.10| ) of the Dirac Hamiltonian 
( |1.16| ) is performed utilizing a spherical box of radius D in order to discretize the mo- 
mentum eigenstates this geometry transfers to the Bethe-Salpeter equation for the meson 
fluctuations. The radial coordinate is then discretized (r^ = Ar{k — 1), k = 1, and 
rj\f = D determines Ar) which transforms the integral equations into matrix equations. 
Then these matrix equations are extended to eigenvalue equations by setting the RHS 
to A(a;)C(r, Co"). The eigenvalue A(ci;) then depends on the frequency u. Finally a; is tuned 
to the eigen-frequency a;, such that X{uJi) = 0. The associated eigenvector represents the 
eigen-wave-function in the descretized form. The bilocal parts of the kernels only involve 
matrix elements of hi ( p.5|) . In these matrix elements upper and lower components of the 
Dirac spinors are coupled because hi is linear in 75. At r = D the unitary transformation 
T (|2.4|) equals unity. Thus the integrands of the matrix elements of hi at r = D are linear 
combinations of terms which are products of an upper and a lower component of the eigen- 
functions of ho at the boundary, ^E'l/dT'l = D). For reasons which will be explained below 
we impose boundary conditions on these eigen-spinors such that the upper components 
always vanish aX r = D It is then obvious that the bilocal parts of the kernels are 
zero as one of the arguments (r or r') lies on the boundary. On the other hand the local 
kernels have finite components a.t r = D. Thus the boundary conditions for the Dirac 
spinors imply boundary conditions for the meson fluctuations as well 



aD,u;) = CAiD,uj) = CBiD,cu) = 0. 



(4.11) 
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For finite D these conditions lead to discretized mesonic modes. 

Since the algebraic expressions for the Bethe-Salpeter equations are by far too com- 
plicated to recognize the appearance of the zero modes we have to verify their existence 
numerically. In order to do so we diagonalize the discretized kernels for cj, = 0. Then we 
compare the wave-functions associated with the lowest eigenvalues with the radial func- 
tions given in eqns ( |4.6| ) and ( |4.7| ). The results are shown in figure 4.1 for the constituent 
quark mass m = 400MeV. The size of the spherical cavity is chosen to be D = 6fm, i.e. 
large compared to the typical extension of the soliton (Ifm). For the rotational zero mode 
we find excellent agreement, while for the translational zero mode a small deviation can be 
observed in the vicinity of the origin. This, however, is due to the fact that the numerical 
computation of the derivative B' near r = is burdened with some small errors. We should 
also mention that the lowest eigenvalue of the Bethe-Salpeter kernel at = is about 
three orders of magnitude smaller than the next to lowest one. I.e. the radial functions 
displayed in figure 4.1 are in fact solutions to the equations (^^ ), ( [4.9|) and ( [4.10|) . Thus we 



have also numerically established the existence of zero modes in the background field of the 
NJL soliton. Furthermore this verification provides an excellent check on our algebraical 
and numerical computations of the kernels which are quite involved. 

Considering figure 4.1 the sizable slope of Cl'™ (^) indicates that the translational zero 
mode has non-negligible overlaps {'H'eT' i''")\'n^^\''^ ^ ^j^"*)) ^i^^ states of large ouj^^. Thus one 
might run into problems concerning the natural cut-off t^^hres = 2m. We will concentrate 
on this when discussing the numerical results for the energy correction in the following 
section. First we have to define the metric for the El and Ml channels. In order to 
take account of the spherical symmetry we define the metric for the Ml channel in the 
discretized coordinate space (rj) via 

Mt'M) = (Arf r^rf /^$r(r,,r;,c.)| (4.12) 



while in the El channel 



(4.13) 

^ " ' ^ U)=U)i 

acts in the 2A^-component vector space spanned by ( V In the discretized coordi- 



C,Bk{^i 

nate space these metric tensors are symmetric matrices A1^[^'^^(u;j). In particular there 

exist transformations l^'s such that the V^s are diagonal. Then the roots ( |2.26| ) are 
defined via the eigenvalue^! \^'^{ui) of M'^^^'^^ 

(v^)^^ (^,) = {v ■ diag yXM). \/Ajv(^.)) ■ V^^^^ (4.14) 
for the Ml and El channels separately. In the latter case the matrices are 2N x 2N 



dimensional. According to ( 2.27 ) the metric is defined into the wave-functions which solve 
the Bethe-Salpeter equation 



uj,) = /\rY,(^/M)^^{uo,)U^,) (4.15) 



*These eigenvalues should not be mixed up with the auxiliary eigenvalues, which were introduced to 
solve the Bethe-Salpeter equation. 
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lEI 
'Ak 
lEI 



aty: 
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AA 

kl 
BA 

kl 



[UJi 



[UJi 



M 
M 



AB 

kl 
BB 

kl 



[UJi 



CAl{uJi 
CBliuJi 



(4.16) 



These modified wave-functions are tlien subject to tlie trivial overall normalization 

l = ArE0r(^.)' and 1 = Ar ^ (^fK^.)' + • 

k k 

Furthermore the overlaps are given by 

{r^^,{r^uj,W'\r,ujf)) = ArY.cl>^\uj{)<t>f''\ujf) 

k 



(4.17) 

(4.18) 
(4.19) 



Here (j)^ ' represent the analogues of in the absence of the soliton. 

Although we are now completely equipped to compute the energy correction /\E we 
postpone this to the next section and rather add some comments on the normalization of 
the zero modes. One can easily show that substituting the rotational zero mode ( |4.2| ) into 
hi (|2.5|) corresponds to 



(4.20) 



with ho being the static Dirac Hamiltonian ( |1.16| ). Since, by definition, the zero mode has 
cUj = one obtains for its normalization from eqn (|2.23|) 



1 



val V- (val|T ■ 5ij|/i)(/i|T ■ <5R|val) 



Cval 



+ 



N, 



^■^(z/|t ■ (5ij|/i)(^|T ■ 8r\u) (e^ 

X 



1/A2 '^^^ A-K Jo 



dx 



|l - sx{l - x) (e^ + euf^ exp -s (^(1 - a;)e^ + a;e^ 
The Feynman parameter integral in this equation can be carried out resulting in 



1 = ^R^Wa, 



where 



a 



^c^vai V- (val|r'^|/i)(/i|rVal) 



ah 



Cval 



(4.21) 



(4.22) 



(4.23) 



^E(Hr1/i)(/i|rV) 



ds 



1/A2 V47rs3 



+ 



is just the moment of inertia for the chiral soliton |^ which actually turns out to be an 
isoscalar, = o?6°'^ . Thus we have shown that the rotational zero mode is normalized 
with respect to the moment of inertia. Analogously one can show that the translational 
zero mode is normalized with respect to the "pushing mass" tensor i^push^"''- This tensor 
is defined by replacing the isospin generators, r° in eqn ( [4. 231 ) by the generators for the 
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infinitesimal translation, ida- In ref.[p^ it has been demonstrated that the pushing mass is 



identical to the energy of the static sohton ( |1.14| ), i.e. -Epush = Ed. Thus the translational 



zero mode is indeed normalized with respect to the classical mass of the soliton. We may 
reverse these results and obtain a possibility to check our metric tensors Ai^^^'^^{r, r', uj = 

0) 

Jd'rJd\'CTir)-M''\r,r',u; = 0)-CTir') = a' dl (4.24) 
J dh j d^r'CTif) ■ M^\r, r\ = 0) ■ Ct"(^') = ^ci 4 (4.25) 



where the matrix structure of hA^^ is not explicitly shown. Eqns ( 4.24| ) and ( 4.25 ) provide 



a possibility to check our calculations in the sense that we first compute the metric tensors 
for cij = and evaluate the integrals on the LHS by substituting eqns ( [4.6|) and ( [4.7| ) 
These results are then compared with the direct evaluations of E^i ( 1.14|) and o? ( [4.23 ) 



Numerically we observe deviations as small of 0.1% when box sizes of the order D = 6fm 
are used and the explicit form in terms of the chiral angle is substituted in the integrals. 
When the solutions to the Bethe-Salpeter equations ai u = are employed to evaluate 
the integrals the error is somewhat larger because these solutions do not exactly match the 
explicit forms at very large radial distances. In any event this error has to be considered 
small and thus provides an excellent verification of the correctness of our algebraical as 
well as numerical manipulations which are rather involved. 

The appearance of the moment of inertia also determined the above mentioned choice 
for the boundary conditions on the Dirac spinors For other boundary conditions, e.g. 
those suggested by Kahana and Ripka ||2^, the moment of inertia is plagued by isospin 
violations of the order 1/D |2^ . 

5. Numerical Results 

In the previous section we have already presented one of our main numerical results: We 
have verified the existence of rotational and translational zero modes. Before concentrating 
on the results for the energy subtraction AE we wish to add a few remarks on the solutions 
in the Ml and El channels when no soliton is present, i.e. G(r) = 0. These solutions are 
important for the overlaps (ry^ ™' |r7*^°)). The Ml channel comes with unit orbital angular 
momentum, / = 1; i.e. a P-wave. Thus the corresponding solutions to the Klein Gordon 
equation which satisfy the boundary conditions ( [4.11| ) read 

Cfrcc(r,^f^) oc ji(g-r) (5.1) 
where the ql make the l^^ spherical Bessel function vanish at the boundary, ji[q\D) = 0. 



Furthermore ul = ym^ + (qj)'^. In the El channel the situation is somewhat more 
involved since S-wave solutions 



a,free(r, uj^'^) cc Uq'^r) and CB,free(r, uf^) = with = ^ml + (g0)2 (5.2) 

as well as D-wave solutions 



CA,free(r,u;f^) oc j2(g«^r) and CB,frce(r, ^) oc ^32{q}r) with uf^ = ^Jml + (qf)^ (5.3) 



exist. We have then computed the solutions to the Bethe-Salpeter equation ( p.l6| ) in the 



absence of the soliton and compared these results with the above suggested solutions to 
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the Klein-Gordon equation. In figure 5.1 we display a typical solution in the Ml channel. 
Here we have chosen = 0. From eqn. ( |2.16| ) we obtain for this solution the eigen- 
frequency lu^^^ = 360.1MeV which reasonably well compares with 363.5MeV as indicated 
by the root q^D of the spherical Bessel function. Except for a small vicinity of r = D 
the radial behavior of our solution to eqn (|2.16|) matches that of the associated spherical 



Bessel function ji{qlr). In figure 5.2 the same comparison is performed for the S- and D 
wave solutions in the El channel. Again the eigen-frequencies which are suggested by the 
solutions of the Klein-Gordon equation are reproduced at the order of 1% and the radial 
dependencies of the solutions to eqn ( |2.16|) reasonably well agree with the corresponding 
Bessel functions. This is especially the case for r < D/2. Since we are interested in the 
overlap with the zero modes, which are well localized (c/. figure 4.1) the deviation from the 
Bessel functions at r D is negligible because this region is not relevant for the overlap. 
Thus we may approximate the solutions to eqn ( p.l6|) in the absence of the soliton by the 
spherical Bessel functions as suggested in eqns ( ^.1|) , (|5^) and (^.3|) . This represents a 
major simplification because solving eqn (|2.16| ) numerically is very time consumingEI. This 



is in particular the case for the El channel since the solutions corresponding to the roots 
g° and qf lie quite close and are difficult to disentangle numerically. 

We thus employ the following procedure to evaluate the matrix elements {ri^'^'\rj^^^): 
We firstly reproduce the zero mode wave-function from the Bethe-Salpeter equation in 
the presence of the soliton and compute the associated metric tensor. This then provides 
the modified wave-functions 0^ ™'(r) in both the Ml and El channels according to eqns 
( [4.15| ) and ( |4.16| ). Next we compute the metric tensors in the absence of the soliton for 
the frequencies ujf'^ determined by the roots of the Bessel functions. Substituting the 
approximations (|5.1|), ( ^.21 ) and ( |5.3| ) into eqns ( [4.15| ) and ( |4.16| ) in turn leads to the 



modified wave-functions (f)^^\r). After normalizing these wave-functions according to eqn 
([4.17]) we are finally enabled to compute the relevant matrix elements as prescribed in eqns 



( [4. 181 ) and ( [4. 19] ). In the course of these calculations we encounter one further problem. 
Due to numerical errors some eigenvalues of the metric tensors turn out to be negative 
yielding eqn ( [4. 141 ) ill-defined. As a matter of fact the absolute values of these negative 



eigenvalues turn out to be about three orders of magnitude smaller than the relevant 
positive ones. In any event we do not expect a numerical accuracy better than about 1%. 
We therefore ignore these negative eigenvalues. This in some sense defines a truncated 
norm. The validity of this treatment can be judged by first normalizing rj subject to 
/ d^r J d?r'ri{r) ■ M.{r, r') ■ r}{r') = 1. Then the modified wave-functions (f> are computed 
with the truncated norm. Numerically we then find the normalization of (p to deviate from 
unity by less than 0.01%. This, of course, justifies the above truncation of the norm which 
turned out to be necessary as a consequence of the numerical inaccuracy. 

We have now completed the presentation of the methods and treatments which serve as 
input to compute the energy correction AE. Needless to mention that we also substitute 



the roots of the Bessel functions for uf'^ in eqn ( 3.23|) . We have already mentioned that, 



a as consequence of the non-confining NJL model, real solutions only exist for uf^ < 2m. 



Due to the treatment in a finite box of radius D the lowest quark energy is wm^ + (tt/D) 



*The computation of the kernels $f | for a given frequency ui takes about 20h-cpu on a HP9000/710 
workstation. In this context it should be remarked that the auxiliary eigenvalues of the Bethe-Salpeter 
kernel for = depend on uj rather strongly. It turns out that uj has to be adjusted with an accuracy 
better than 0.1% in order to solve the Bethe-Salpeter equation. 
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Table 5.1: The quantum corrections to the sohton mass due to the rotational zero mode. 
The size of the spherical cavity is D = 6fm. 



m(MeV) 


= 

400 500 600 


= 135MeV 
400 500 600 


s 

AE(McV) 


0.88 0.94 0.95 
-201 -274 -290 


0.83 0.89 0.91 

-244 -297 -323 



which increases the threshold energy to cu^^ = 2^m? -f {n/Dy in the Ml channel while 

in the El channefl = \Jm^ -f [n / DY + s^rn? -f (^i)^. We therefore truncate the sum 
( p.23|) accordingly. In order to judge this truncation we define the sum of overlaps for the 
zero modes 

5= E |(^7^-"^-(r')|^(°)(r',^f))f (5.4) 

J — 

which should approach unity if the model were insensible to the truncation. Note that the 
number of meson modes which lie below cu*^ decreases as m^r increases. 

In table ^]T] the results for the rotational zero mode are presented. We see that the 
sum of overlaps S is about 0.9 for all sets of parameter used. This number, of course, 
increases with the constituent quark mass m because the threshold grows proportionally. 
Furthermore S is sufficiently close to unity in order to conclude that the energy subtraction 
/XE ^ -(250 - 300)MeV is reliable. Actually this is about lOOMeV smaller than the 
corresponding value in the Skyrme model |jl^. We also observe that the mass correction 
due to the rotational zero mode is about (30-40) MeV lower in the chiral limit than for the 
physical value of the pion mass, m.^ = 135MeV. 

Table |5.2| contains the results for the translational zero mode. We see that the sum of 
overlaps is significantly smaller than for the rotational zero mode (table p?TD . Here it hardly 
reaches 0.5. This is due to the fact that the slope of the profile function corresponding to the 
translational zero mode is enhanced compared to that of the rotational zero mode. Thus 
a sizable number of Fourier components with non-vanishing overlap lie beyond the quark- 
antiquark threshold. Hence we have to interpret the resulting AE ^ —(100 — 200)MeV 
as a lower bound for the energy subtraction originating from the translational zero mode. 
Nevertheless we can extract some qualitative statements from the results listed in table 
|5.2|. First of all we observe that the D-wave contributions are dominating. This result is 



also found in the Skyrme model |1J]. Furthermore we see that the value of the pion mass 
has only little infiuence on the S-wave contribution while the absolute value of the D-wave 
contribution decreases with increasing pion mass. This is in contrast to the rotational zero 
mode. Thus the net effect of varying the pion mass is somewhat mitigated. 

Next we wish to estimate the contributions of the scattering states, Ui > itlt^, to the 



energy correction Ai? as described by eqn (|3.23| ). In order to do so we first have to 



construct the corresponding solutions to the Bethe-Salpeter equation ( p.l6|) . In table ^ 



the frequencies of the lowest solutions are compared to their analogues in the absence of 
the soliton. In the El channel we distinguish between solutions which are dominantly S- 



^^Note that in the El channel hi has vanishing matrix elements in the grand spin zero subsystem of the 
quark modes. 
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Table 5.2: The quantum corrections to the sohton mass due to the translational zero mode. 
The contributions stemming from the S{1 = 0)- and D{1 = 2)-waves are disentangled. The 
size of the spherical cavity is D = 6fm. 



m(MeV) 


= 

400 500 600 


= 135McV 
400 500 600 


s 

AE/=o(MeV) 
AE/=2(MeV) 
AE{MeY) 


0.41 0.41 0.61 

-18 -22 -28 
-127 -140 -207 
-145 -162 -235 


0.28 0.37 0.46 
-12 -22 -28 
-82 -128 -187 
-94 -150 -215 



Table 5.3: The contribution of the first scattering states to the quantum corrections of 
the soliton mass in the channel of the rotational (Ml) and translational {El) zero modes. 
The scattering state under consideration is labeled r]. The entry "max(..)" represents the 
contribution to S by the dominant term. The constituent quark and the pion masses are 
600MeV and 135MeV, respectively. The size of the spherical cavity is D = 6fm. 





Ml 


El 


/ = 1 


I = 


/ = 2 


u;(°)(MeV) 


202 


171 


235 


u;(MeV) 


210 


176 


240 


max ((ryf^lry)^) 


0.96 


1.00 


0.99 


S 


0.99 


1.00 


1.00 


AE(MeV) 


-1.9 


-0.6 


-0.5 



or D-waves. We establish that the scattering states always lie slightly above the associated 
states for B = 0. We furthermore observe from table |57^ that the sum over the overlaps S is 
strongly dominated by only one term. As a matter of fact it is always the one which happens 
to have the same number of knots in the radial part of the wave-function associated with 
scattering state under consideration. This dominance causes S to very closely approach 
unity. Thus we conclude that the contribution of the low-lying scattering states to AE 
of the order of only a few MeV is very reliable. We have found the same result for the 
second scattering state in the Ml channel. For the El channel the extraction of higher 
scattering states is somewhat troublesome because the states which are dominantly S- or 
D-waves are almost degenerate. Nevertheless the results found so far for the scattering 
states suggest that their contribution to AE is almost negligible. The fact that these 
contributions turn out to be negative is a consequence of the expansion (|3.21|) when the 



"perturbation" V is attractive [^. These results qualitatively agrees with those obtained 



in the Skyrme model [14 



We are finally enabled to present our predictions for the masses of baryons in the NJL 
soliton approach. For the projection onto good quantum numbers of spin and isospin we 
employ the well-established semi-classical cranking approach ||^. This yields the mass 
formula for a baryon of spin J 

M = E,, + AE+^^^. (5.5) 
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Table 5.4: The predictions for the masses of the nucleon (A^) and A-resonance. The 
empirical data are 939MeV and 1232MeV, respectively. 







m^r = 






= 135MeV 


m(MeV) 


400 


500 


600 


400 


500 


600 


Eci(MeV) 


1212 


1193 


1166 


1250 


1221 


1193 


Aii;(MeV) 


-346 


-436 


-525 


-338 


-448 


-538 


a2(i/GeV) 


6.26 


4.73 


3.87 


5.80 


4.17 


3.43 


M7v(MeV) 


926 


836 


738 


976 


863 


764 


MA(MeV) 


1166 


1153 


1126 


1236 


1223 


1201 



According to the above discussions AE is understood as the sum of the zero mode con- 
tributions listed in tables |5.1| and f).2[ In eqn (|5.5|) refers to the moment of inertia 



( ^ , ^) [0]. The individual pieces in eqn are of the orders 0{Nc), C(l) and 

0{1/Nc), respectively. Here we have ignored the quantum corrections at order 0{1/Nc)- 
These have been estimated in the Skyrme model to be less than lOOMeV [M . Hence their 



contribution seems to be less than the uncertainties in our estimates for AE. 

We indeed find that the numerical results roughly follow this 1/Nc counting pattern. 
Unfortunately the results shown in table ^.4| suggest that the overall prediction for the 
baryon masses in the NJL soliton model underestimates the empirical value for the nucleon 
mass (939MeV). Only for the constituent quark mass m ^ 400MeV a good description is 
obtained. In that case, however, S is as small as 0.4 in the El channel (c/. table ^.2|) . 
Then, of course, the question of reliability cannot be answered unambiguously. Although 
one may assume the point of that this is a shortcoming of the NJL model in general rather 
than only for the meson fluctuations. Hence one would have to consider the results listed in 
tables ^TT] and |5.2| seriously and conclude that the masses of the nucleon and A-resonance 
are reasonably well reproduced for m ~ 400MeV. As m is further increased the classical 
soliton energy E^i as well as the quantum corrections decrease; thereby underestimating 
the nucleon mass. Simultaneously the 1/Nc counting pattern breaks down since also the 
moment of inertia, a^, decreases leading to a A-nucleon mass splitting of the order of 
AE. These two quantities are, however, supposed to differ by one order in Nq. In turn 
this renders the prediction for the mass of the A-resonance almost independent of the 
parameters and in reasonable agreement with experimental value of 1232MeV when the 
physical value of the pion mass is adopted. 

6. Conclusions 



In the present paper we have studied the quantum corrections, AE, to the classical 
mass. Eel, of the chiral soliton in the NJL model. In a first step we have investigated 
the formal structure of the Bethe-Salpeter equation for pionic fluctuations in the soliton 
background. As the soliton is static and the eigen-modes of this equation appear in pairs 
±uji we were able to demonstrate that the energy operator in the Fock space of the meson 
fluctuations is identical to that of a harmonic oscillator; although the time derivative 
operator appears at all (even) orders in the coordinate space representation of the Bethe- 
Salpeter equation. This result has allowed us to adopt the expression for the quantum 



correction which was previously obtained in the Skyrme model |14|. According to our 



examinations this expression is applicable to all static soliton models in 3 -|- 1 dimensions as 
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long as the eigen-modes occur in pairs. However, here the overlap matrix element between 
solutions of the Bethe-Salpeter equation with and without the soliton present, has turned 
out to be more involved than in the Skyrme model in particular because the solutions in 
the absence of the soliton do not exactly obey the Klein-Gordon equation. Furthermore 
the sum over the free eigen-modes had to be truncated since these modes are unstable 
against the decay into a quark-antiquark pair once the corresponding eigen-frequency has 
exceeded the threshold, 2m. 

Two types of eigen-modes exist in the presence of the soliton: The bound zero modes 
and the scattering states. We have made plausible that the latter type only contributes 
very little to /\E. Hence the scattering modes have been discarded and /\E has completely 
been approximated by the zero mode contribution. This procedure is also justified by the 
similar calculations in the Skyrme model [p!4 



Before actually computing the quantum corrections we have, for the first time, verified 
the existence of rotational and translation zero modes in the background of the NJL soliton 
by explicit construction. For the rotational zero mode the truncation caused by the non- 
confining nature of the NJL model of the mode sum does not represent a serious problem 
because the sum of overlaps approximates unity fairly well. This in turn leads to a reliable 
estimate for the quantum correction due to the rotational zero mode of about — 250MeV. 
In the case of the translational zero mode the situation is somewhat worse because the 
overlaps sum up to only about 0.5. In this sense one has to regard /\E ^ — ISOMeV only 
as a bound. 

One might, however, adopt a different point of view and consider the above quoted data 
as the actual results. In any event the NJL model is not well defined for frequencies above 
the quark-antiquark threshold. All other quantities {e.g. the classical mass Ed) might 
undergo significant changes as well once the model is improved to avoid this problem. In 
this interpretation we have observed that the NJL model predicts the masses of the nucleon 
and the A resonance fairly well for a constituent quark mass m ~ 400MeV. Even the I/Nq 
counting scheme Ed ~ NqAE ~ N^^Mi^ — Mn) seems to operate then. Upon increasing 
the constituent quark mass this is no longer the case and the estimate for the nucleon mass 
turns out to be somewhat too small while the prediction for Ma remains almost unaltered. 
The smallness of the absolute value for the nucleon mass appears to be connected to the 
fact that Ed decreases as the constituent quark mass increases. This seems to be special 
to the NJL model without the isoscalar-vector cj-meson included. Recently it has been 
demonstrated that the proper incorporation of this field indeed yields a classical energy 



which increases with m Thus one might suspect that the associated extension of the 
NJL soliton model leads to an even more reasonable description of the nucleon mass for 
larger constituent quark masses. This would in a sense be more reliable because the sum 
of overlap matrix elements would come closer to unity. 

As a side-product we have been able to define a metric for the overlaps of the meson 
states with and without the soliton present. As has previously been pointed out in the 
context of the Skyrme model |31| this metric plays an important role for the computation 



of the momentum dependent pion-nucleon form factor gnNNiQ"^)- In order to compute 
QirNNiff) one nai'ely would Fourier-transform sinB(r). Commonly in soliton models this 
procedure underestimates the cut-off Attatat ~ 1.6GeV, which is defined via 



9nNN{q^ < 0) a: 



2 ^2 



ttNN 



(6.1) 



by a factor two or even more [^]. This may be improved by the proper incorporation of 
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the metric tensor because the Fourier-transformation of sin0(r) basically represents the 
projection of the rotational zero mode onto free pion states. Investigations in this direction 
in the context of the NJL model are subject to future studies. 
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Figure captions 



Figure 4.1 

Comparison of the numerical solution of the Bethe-Salpeter equations (4.8,4.9,4.10) with 
the analytic form for the zero modes as given in eqns (4.6) and (4.7). The lattice points 
are indicated. 

Figure 5.1 

Comparison of the numerical solution of the Bethe-Salpeter equations (2.16) for = 
with the corresponding spherical Bessel function in the Ml channel. Both curves refer to 
the wave-function of the second excited state. The normalization of the radial functions 
is arbitrarily chosen. The lattice points are indicated. 

Figure 5.2 

Comparison of the numerical solution of the Bethe-Salpeter equations (2.16) for G = 
with the corresponding spherical Bessel function in the El channel. Left: The S'-wave 
solution with Cb(^) = 0- The normalization of the radial functions is arbitrarily chosen. 
The corresponding solution of the Klein-Gordon eqn has lo^^^ — 314.2MeV. Right: The 
D-wave solution. Here the linear combination Ca('") — Cb('")/3 vanishes. The Klein-Gordon 
eqn yields the eigen-frequency 303.2MeV. 
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